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Abstract 

We compute the contribution from the twist-3 fragmentation function for light hadron production 
in collisions between transversely and longitudinally polarized protons, i.e., pip —>■ hX, which can 
cause a double-spin asymmetry (DSA) Alt- This is a naive T-even twist-3 observable that we 
analyze in collinear factorization using both Feynman gauge and lightcone gauge as well as give 
a general proof of color gauge invariance. So far only twist-3 effects in the transversely polarized 
proton have been studied for Alt in P^P hX. However, there are indications that the naive 
T-odd transverse single-spin asymmetry (SSA) A^ in pip h X is dominated not by such distri¬ 
bution effects but rather by a fragmentation mechanism. Therefore, one may expect similarly that 
the fragmentation contribution is important for Alt- Given possible plans at RHIC to measure 
this observable, it is timely to provide a calculation of this term. 

Keywords: twist-3 effects; spin asymmetries; fragmentation 


1 Introduction 

Spin asymmetries in various hard processes have brought a new perspective to high-energy perturbative 
QCD theory and phenomenology. Transverse single-spin asymmetries (SSAs) Ayr, which are naive T- 
odd observables, were first explored in the mid-1970s in pip —vr A at Argonne National Lab [1] and 
pBe —7> A^A at FermiLab [2]. Both of these measurements led to strikingly large effects that were 
unexplainable in the naive parton model [3]. Experiments continued at FermiLab in the 1990s for 
pip TT X ^ and most recently at AGS [5l[6] and RHIC [THH] for pip —> {vr, K, jet}X. All of 
their measurements likewise produced substantial transverse SSAs. On the theoretical side, it was 
realized in the 1980s by Efremov and Teryaev that if one went beyond the simple parton model 
and included (collinear twist-3) quark-gluon-quark correlations in the nucleon, then there was the 
potential to generate these large effects [15]. A systematic approach was then developed by Qiu and 
Sterman in the 1990s that presented the collinear twist-3 factorization framework |161tl8] with the 
expectation that one would be able describe transverse SSAs within this perturbative approach. Later 
a solid foundation was given to this formalism in |191I20] that proved the cancelation among gauge- 
noninvariant terms and led to an expression for the twist-3 cross section in terms of the complete set 
of the twist-3 quark-gluon-quark correlation functions. Over the last decade, several other analyses, 
including those for the extention to twist-3 fragmentation functions [2IH24] and three-gluon correlation 
functions |25] . furthered the progress of this formalism — see also |261l3T] and references therein. 
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For many years the main assumption was that these transverse SSAs were due to effects inside 
the transversely polarized proton, in particular those embodied by the so-called Qiu-Sterman function 
Tp However, a fit of the QS function to A]\f data led to a result that was inconsistent with 

an extraction of the Sivers function [32] from SIDIS, which has a model-independent relation to 
Tp |33|, and became known as the “sign mismatch” crisis [M]. An attempt to resolve this issue through 
more flexible parameterizations of the Sivers function proved unsuccessful [35], and, by looking at Atv 
data on the target transverse SSA in inclusive DIS [361IF7] . it was argued in fact that the QS function 
could not be the main cause of Aw [38] . This led to a recent work that examined the impact of 
fragmentation effects from the outgoing hadron [32] based on the analytical calculation in Ref. [23j . It 
was determined that this fragmentation term could be the dominant source of Aw in p^p —>■ irX |39j . 

In addition to Aw, there is another twist-3 observable in proton-proton collisions that can give 
insight into quark-gluon-quark correlations in the incoming protons and/or outgoing hadron. This is 
the longitudinal-transverse double-spin asymmetry (DSA) App, which, unlike Aw, is a naive T-even 
process. The classic reaction for which this effect has been analyzed is App in inclusive DIS (see m 
for recent experimental results on this observable). This asymmetry has also been studied in the Drell- 
Yan process involving two incoming polarized hadrons |41H44j : in inclusive lepton production from 
VF-boson decay in proton-proton scattering [45]; for jet production [46] and pion production m in 
lepton-nucleon collisions; and for direct photon production [IS], jet/pion production [l9], and D-meson 
production [50] in proton-proton collisions. 

Of these works on App, only in Ref. [U] for ip^ —>• ir X was the twist-3 fragmentation piece 
calculated (we will see the structure of that result persists in our computation), whereas the frag¬ 
mentation term for p^p —>• vr A has never been studied. Like with Aw, there is no reason a priori 
that this piece cannot be important or perhaps dominant in the asymmetry. Given the possible plans 
by the PHENIX Collaboration at RHIC to measure App for pions [51^1 . we feel a calculation of the 
fragmentation term for this final state is needed at this time. Furthermore, like prior research in the 
literature [2TIl24l l47 l l54 [ l55j . this work will continue to establish/verify the theoretical techniques used 
in collinear twist-3 fragmentation calculations. The remainder of the paper is organized as follows: in 
Sec. Owe introduce the twist-3 fragmentation functions relevant for spin-0 hadron production. Next, 
in Sec. [3| we discuss the calculation of the polarized cross section formula for App. Finally, in Sec. O 
we summarize our work and give an outlook. A general proof of the color gauge invariance of our 
result is given in the Appendix. 

2 Twist-3 fragmentation functions for spin-0 hadrons 

We now define the set of twist-3 fragmentation functions relevant for spin-0 hadron production. The 
quark-quark matrix element gives two purely real twist-3 functions, which read 

/ ^e-^Ho\rim\Ph;X){Pp;X\^]{Xwm 

X 

= Zp(z) + tpiz) + ■■■, ( 1 ) 

where ijji is a quark field with spinor index i, and we use the simplified notation 

(with eoi 23 = +!)• The color indices are summed over and divided by the number of colors N = 3. 
The scale is used to make the functions dimensionless and is on the order of the nucleon massU 

^We mention that a clear Alt asymmetry has already been seen by the Hall A Collaboration at Jefferson Lab in 
SIDIS [52] and IrA ^ n X 

'^Mx is the scale of nonperturbative chiral-symmetry breaking (CSB), which is said to be on the order of the nucleon 
mass IGeV). We have chosen this scale instead of the light hadron mass Mh since twist-3 functions representing 
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The vector w'^ is light-like (u)^ = 0) and satisfies Ph-w = 1. We will suppress the gauge-link operators 
throughout for simplicity. 

Next, we introduce the so-called F-type quark-gluon-quark twist-3 fragmentation functions. We 
can dehne two independent functions as 


^ ^i)(0|V'f(0)|P,;X)(P,;X|^J(Au;)5F“-(/iu;)|0) 


( 2 ) 




(3) 


where is the gluon field strength tensor. We note that both Ep{zi,z) and Ep{zi,z) in 

general are complex functions. The correlator Ep{zi,z) has support on 1 > z > 0 and > z, while 
Ep{zi, z) has support on ^ ^ > 1, ^ < 0, and ^ > 0 [MllSG] . 

We can consider the so-called F-type twist-3 fragmentation functions Ep{zi,z) by replacing 
gF°‘'^{fj.'w) in ([2]) with a covariant derivative D'^{g,w) = 9“ — igA°^{gw). However, Ep{zi,z) and 
Ep{zi, z) can be related through the identity, 

where e{z) is another twist-3 fragmentation function that is pure imaginary and defined as 


1 

N 




{0\[oow,0]'il;'^{0)\Ph-,X){Ph;X\'il;'^-{Xw)[Xw,oow]\0)^° 


M^c 

2z 




(5) 


Note that we have restored the gauge links [a, b] in order to emphasize that d°‘ acts on the A dependence 
in both the quark field '4)'j{Xw) and the gauge link [Arc, oorc]. The F-type function Ep{zi,z) has 
another relation associated with the QCD equation of motion. 


dz 

P T 


dzi ^h/qf 


Ed'‘{zi,z) = ef {z) + i e'' (z). 


h/q. 


( 6 ) 


By combining Eqs. dH, ([6]) we can eliminate the F-type function and obtain 


T“^P 

Jz zf 


1 


Xjzi - Ijz 


) Ey\z^,z) + zP^I<i{z) = e'll\z) Piell\z) . 


(7) 


The real and imaginary parts of the above relation respectively give 





dzi 


P 


l/zi - Ijz 


E 


h/q,'St 


(zi,z) = ef'^®(2;) 


( 8 ) 


helicity flip effects are due to nonperturbative CSB whereas Mu for the pseudoscalar mesons represents the explicit CSB 
due to the quark mass. 
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(9) 


■i: 


dzi 


P 


1/zi -1/z 


^ Ep^‘^’''{zi,z) + ze’^^‘^’'^{z) = e'l^'^{z) 


where 3? (9) indicates the real (imaginary) part of the function. It was shown that Eq. Q ensures the 
gauge invariance of the polarized cross section formula in the case of the transverse SSA in SIDIS [24] . 
We will show that Eq. ([8|) plays the same role in the case of the longitudinal-transverse DSA in 
proton-proton collisions. 


3 Calculation of the polarized cross section for Axt 

We consider the polarized cross section for the production of a light (spin-0) hadron from the collision 
between a transversely polarized proton and a longitudinally polarized proton, 

p(P,Sx)+p(P',A)^/i(PO + ^, (10) 

where the momenta and polarizations of the particles are given. The first non-vanishing contribution 
to the cross section reads 


daiPhU S±,A-) = H ^ fa/A{3) ® fb/B{2) ® DcIc{2) 

+ H' iSl fa/A{2) ^ fb/B{3) ^ -^C/c(2) 

+ H" (g) fa/A{2) ^ fb/B{2) ® -^C/c(3) ; (H) 

where a sum over partonic channels and parton flavors in each channel is understood. The labels 
on the functions indicate the parton/proton (or hadron/parton) species and the twist (e.g., fa/A(3) 
denotes a twist-3 correlator associated with parton a in proton A). These functions are convoluted 
with hard factors H, H', H", which are different for each term. The first term in m was already 
calculated in Ref. [l9]. With regards to the second term, for the case of the transverse SSA (where 
B is now unpolarized), which involves chiral-odd twist-3 unpolarized distributions, the authors of 
Ref. |57j demonstrated this part is negligible because of the smallness of the hard scattering coefficients. 
However, in Ref. |43j the authors found in Drell-Yan for Alt this second term, which involves the 
chiral-odd twist-3 longitudinally polarized distribution, can be as large as the first (chiral-even) term. 
Therefore, this second term should be analyzed in the future in order to have a complete result for this 
observable. We will now compute the third term, which involves the chiral-odd twist-3 fragmentation 
functions introduced in Sec. [2] coupled to the (chiral-odd) transversity function hi{x) and the helicity 
distribution gi{x), both defined in the standard way 

I S\fPlmUXn)\P, S) = i fiix) + MAPh stix) + h\{x)\ , (12) 


f Ap^lp,S\F“"{0}F^''{\n)\P,S) 
J 27r 



a'f fi{x) + gl (x) 


(13) 


where the unpolarized distribution fi{x) has been included for completeness. The vector is light-like 
and satisfies P ■ n = 1. 

The techniques for calculating the complete fragmentation term in the collinear twist-3 framework 
have been laid out in Refs. |23ll24l[55] . In particular, the work in |23] can be used for deriving the 
result in lightcone gauge, whereas that in [24] can be employed if one chooses Feynman gauge. For 
the former, one can make straightforward changes to the calculation in Ref. [23] in order to obtain the 
result for p^p -A- hX (since a similar process {p^p -A hX) was computed there). The channels that 
one must consider are qg -A qg, qq' -A qq', qq -A qq, qq -A qq, qq' -A qq', and qq -A qq (see Figs. 
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Figure 1: Feynman diagrams for the qg —>■ qg channel. The circled cross indicates the parton that 
fragments. The 9 graphs (ignoring the dots) lead to the hard factor for ei{z) while those with the dots, 
which represent coherent gluon attachments from the parton line to the fragmentation correlator, give 
the hard part for Ep(zi,z). The Hermitian conjugate (H.c.) graphs are also taken into account. 


plus all the antiquark fragmentation channels found through charge conjugating the aforementioned 
ones. We note that because of the T-even nature of this observable, the structure of our cross section 
takes on a different form than Ai^. For example, one no longer has contributions from ei{z) or e{z)-, 
instead, the quark-quark piece only involves ei{z). Likewise, for the quark-gluon-quark part, one 
receives contributions from Ep{zi,z) instead of Ep{zi,z). The hard factors for Ef{zi,z) (see Fig. SD 
vanish, as was also found in Ref. [23] with Furthermore, unlike Aw, it turns out that the quark- 
gluon-quark hard factors are independent of zi. Thus, one can use Eq. (I8|) to write this piece in terms 
of 61 ( 2 ;), and the entire cross section then only involves this function. Such a simplification was also 
noticed in Ref. [17| for £p^ tt X. We then find the fragmentation term in the polarized cross section 
relevant for Alt to be 


(PPh 




S 


^ i a,b,c 

The hard factors at are given by 


f ^g\{x') f — h1{x) ai5{s+t+u). (14) 

A ^ Jo X Jo X 


<7, 




^s'^ + u^ 11' 

i'^u u 


'Xqq’-i-qq' — I 1 


iv2 1 ’ 


^qq^qq — 1 1 


1 


s I 1 s{s — 2i) 


^ nn—vnn — I 1 


, ^qq^qq 


1 


s 1 u — 2t 
¥ ^ N 2P 


(15) 


^qq^qq-j^[^ iV2 j 2u ' 
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Figure 2: Feynman diagrams for the qq —>• qq channel. The first graph gives the qq' —)• qq' channel. 





Figure 3: Feynman diagrams for the qq —>■ qq channel. The first graph gives the qq' —>• qq' channel. 
Note that the qq —?> qq channel is found by interchanging the two incoming parton lines. 

The Mandelstam variables for the process are defined as S = {P + P')'^, T = {P — Ph)'^, and U = 
{P' — Ph)'^, which on the partonic level give s = xx'S, t = xT/z, and u = x'U/z. For the antiquark 
fragmentation channels we find (cf. [23]) = ^ab^cd, where &ab^cd are given in (fT5]) . We also 

calculated the polarized cross section in Feynman gauge using the procedure of Ref. |24] and found 
agreement with Eqs. (I14p . (I15p . A general proof of color gauge invariance can also be found in the 
Appendix. 

4 Summary and outlook 

Transverse SSAs A^r in single-inclusive processes (e.g., p^p —>■ hX) are twist-3 observables that have 
been an intense topic of research for close to 40 years. Large effects have been found that still have an 
unclear origin. Recently it has been shown that the fragmentation term in collinear twist-3 factoriza¬ 
tion could be the main cause of Ajv in p^p tt X at RHIC |39] . In addition, another twist-3 reaction 
exists that can also lead to information on quark-gluon-quark correlations in protons/hadrons: the 
longitudinal-transverse DSA Alt in P^P hX. Already two related observables, Alt in SIDIS [52] 
and in £n^ —> ttX [53], have been measured and nonzero effects have been found. However, RHIC, 
with the only source of (independently manipulated) polarized proton beams in the world, has yet to 
explore Alt in p^p hX despite measuring asymmetries for every other combination of proton spins. 
Just recently, though, the PHENIX Collaboration has put forth plans to make this measurement m- 
Some work has been done previously on Alt that looked at the twist-3 effects in the polarized pro¬ 
ton [49]. Motivated by the potential of twist-3 fragmentation effects to dominate Ajsr [39], and with no 
reason a priori that they should be small, we have computed this term now for Alt- We found that 
the entire result can be written in terms of a single (twist-3) quark-quark fragmentation function e\{z) 
and confirmed this using two different gauges. In the future we plan to perform a detailed numerical 
study of Alt in P^P —>■ tt A in order to further encourage an experiment at RHIC. 
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Figure 4: Feynman diagrams for qg and qq induced channels where both quarks entering the frag¬ 
mentation are on the same side of the cut. These lead to the hard factors for Ef{zi, z). Such graphs 
cancel after one sums all contributions. 
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Appendix: General proof of color gauge invariance 

We want to show that our calculation of the qg qg channel satisfies the Ward-Takahashi identity 
(WTI): 

{xP + x'P' - Ph/z)P = 0 , (16) 

where is the partonic hard factor with the polarization tensor for the hnal unobserved gluon 

dP^{xP + x'P' — Ph!z) removed. The momenta and Lorentz indices are shown in Fig. 5(a). In our 
computation, we wrote da’^P^^P as follows (cf. Fig. [T]): 

, (17) 
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Figure 5: Feynman diagrams that enter our general proof of color gauge invariance: (a) quark-quark 
and (b), (c) quark-gluon-quark graphs. Note that our calculation does not include in (b) the subset 
of diagrams given in (c), but rather these topologies are taken care of by those in (a). See the text for 
details. 


where the superscripts correspond to the diagrams in Fig. [5] (L indicates the explicitly shown graph 
and R its H.c.). One sees immediately that 

{xP + x'P' - Ph/zY = 0 , (18) 

since from the WTI for qg —>■ qg scattering we know 

{xP + x'P' - Px/zf {>t'‘, } = 0 (19) 

if all particles attaching to the respective blob are external (on-shell) lines. The Dirac projections 
used for the correlators in dY^^^ allow these connecting particles to meet this criteria]! 

Therefore, we now have 

(^P + YP' - Ph/zY da^r^^ = [xP + x'P' - Ph/zY (dafj - da^Y^) 


~ {xP + x'P' - Ph/zY 
'Tr 


X i - I ’ 


(A/^)757r 


{PhIz) 


2" PM 


-I- Tr 




( 20 ) 


One can show that 

[xP + x'P' - PhjzY MY. Y 


®The diagram in Fig. 5(a) does not satisfy this requirement because the fragmentation correlator is projected out 
with 1 . 










AluPn 


[fh/z] 


i{xf^x't) , PhiJilz ,Mcrj.d 


■1^T^T^ + 


jodc rpd 


= [f>h/z) 

where are the structure constants and T°' the generators of SU{3). This leads to 

+ x>p' _ p^/^)P ^ _1 pr'P,^ {^hlz)75lrM' 


( 21 ) 


+ Tr 




( 22 ) 


We can simplify this expression further through use of the identities 


' 15 (^ 111 ' 


1 

2 







(23) 


and obtain 


{xP + x'P' -Ph/z)Pdall-^‘i<^ 


1 


■Tr 


2 

+ Tr 

= 0 . 


$xfl^Ml^^{fhlz){fhi^ - 




'ucb 


(24) 


Thus, the qg —> qg channel satisfies the WTI. 
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